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Manipulation of a single magnetic atom using polarized single electron transport in a 

double quantum dot 
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We consider theoretically a magnetic impurity spin driven by polarized electrons tunneling through 
a double quantum dot system. Spin blockade effect and spin conservation in the system make the 
magnetic impurity sufficiently interact with each transferring electron. As a results, a single collected 
electron carries information about spin change of the magnetic impurity. The scheme may develop 
all electrical manipulation of magnetic atoms by means of single electrons, which is significant for 
the implementation of scalable logical gates in information processing systems. 
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I. Introduction 

Magnetic atoms are critical spin systems which have 
potential applications in data storage and quantum in¬ 
formation processing ^ -3 In particular, using electrons to 
manipulate magnetic atoms is a natural step towards the 
implementation of scalable memory units for future inte¬ 
grated circuits. In dilute II-VI semiconductor quantum 
dots (QD), interaction between a manganese (Mn) atom 
and a carrier can be effectively described with the sp — d 
exchange interaction^ - — Based on the impurity-carrier 
coupling, electrical control of the single magnetic atom 
is feasible by injecting different charges in the magnetic 
atom doped QD— 

As shown both in experiment— and theory/ electrons 
can directly tunnel through a Mn atom by taking its spin 
states. The ground state and excited states of the Mn 
atom can be identified with the current, since they sup¬ 
port different conductance. However, in the tunneling 
from a STM tip to Cu^N surface through an individual 
Mn atom, the Mn spin spontaneous relaxation is more 
frequent than excitation by tunneling electrons— Com¬ 
pared with the case in high dimensional bulk material, 
the life time of the Mn spin is much longer in a QD— 
Therefore, coherent electrical manipulation of the Mn 
spin is possible in low dimensional nanostructure. In a 
QD doped with a single Mn atom, charge and conduc¬ 
tance of the single electron tunneling can be related to 
the spin state of the Mn atom— Whereas, it is hard to 
exactly connect a quantum state of the magnetic atom to 
a single electron in the above QD system, which should 
be a problem required to be solved for future quantum 
information processing. 

In this paper, we propose a scheme for all electrical 
manipulation of the magnetic impurity spin, scaling the 
number of driving electrons down to one. To this end, we 
consider two inter-coupled semiconductor QDs, in which 
one QD contains a single magnetic impurity and couples 
to a spin polarized electron source. The other QD is lo¬ 
calized in a homogeneous magnetic field and connected 
to a normal conductance, playing a role of a spin filter. 
In previous study/^— — electrons transport through the 
magnetic atom without any qualification to their spins, 
as a result, interaction between the magnetic atom and 


electrons are very weak. In contrast, the QD spin filter 
induces spin dependent tunneling, which makes sure each 
electron would be completely flipped by the magnetic 
atom before it passes through the filter. Therefore, each 
collected electron can be correlated with the change of 
spin state in the magnetic atom. There are several facts 
that are very beneficial for the realization of our scheme. 
First, coupled two QDs can be fabricated, doping a single 
magnetic ion in one of them— Second, spin life time of 
an electron in II-VI semiconductor quantum dots can be 
long enough. Electron relaxation time in a similar system 
was reported to be 50 ns in a previous work * 11 In a QD 
imbedded a magnetic atom, a longer life time of around 
1 p,s was predicted : 12 Third, Magnetic atom such as Mn 
impurity with relaxation time from 1 ps to 0.4 ms was 
observed in experiment —?’ 13 The time is longer than driv¬ 
ing time of the magnetic atom, which will be shown later 
in the present work. Forth, the technique for real time 
detection of single electron tunneling has been well de¬ 
veloped recently— - — Fifth, polarized electron current is 
available from several sources, for instance, QD spin split¬ 
ter under the premise of local magnetic field , 18 ’ 19 ferro¬ 
magnetic leads , 20 ’ 21 graphene or carbon atom wires : 22 ’ 23 
Recently, nuclear spins of donor atoms such as phospho¬ 
rus and 29 Si have been coherently controlled and read 
using bounded electrons in these donors and the spin to 
charge conversion. Electromagnetic field was applied to 
initialize these nuclear spins. 2 - - — Comparing with these 
early works, the main advance in our present protocol is 
that the magnetic atom would be manipulated all elec¬ 
trically. In other words, a magnetic atom can be con¬ 
trolled only using a single electron transistor without the 
application of any electromagnetic field. Even the initial¬ 
ization of the impurity spin can be progressed using the 
single electron tunneling in principle. 

II. MODEL AND THEORY 

Our model is illustrated in Fig. [0 The two QDs are 
denoted by dot 1 and dot 2 with ground orbital levels 
£1 and £ 2 , respectively. Both electron polarization in the 
left lead and the external magnetic field B that applied 
on the dot 2 are assumed to be parallel to the QD growth 
direction 3 . Electrons injected from the left lead into the 
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FIG. 1: (Color on line) Schematic illustration of the principle 
in our model. In (a) - (d) the applied magnetic field is in 
down direction and up polarized electrons are injected into the 
double QD. In (e)-(h), the external magnetic field is changed 
to be up direction and down polarized electrons are injected. 


dot 1 are coupled to the magnetic atom by the ferro¬ 
magnetic Heisenberg type spin exchange interaction. We 
describe spin of the magnetic atom with mean value of 
the spin along z direction (M z ), where M is the mag¬ 
netic atom spin operator and the bracket indicates aver¬ 
age over quantum states of the system. Bias voltage and 
the magnetic field is tuned that just the lowest levels of 
the two QDs fall within the bias window g k > £i, (or 
e t)> /Hi, where the indexes j., f indicate electron states 
with spin up and down, respectively. In addition, the 
intra-dot Coulomb blockade energies Ui, U 2 correspond¬ 
ing to the dot 1 and the dot 2 are assumed to be much 
larger than other energy scales, which yields only single 
electron occupation is involved in either of the QDs. It 
remarkably simplifies our model and calculation. 

In Figs, QJa)-(d), the magnetic field in this config¬ 
uration is applied along the — z direction with value 
B = (0,0, —B). It means, in dot 2, the ground state 
level is £ 4 , = £2 — g* ^bB/2 and the first excited level is 
£ 4 - = £2 + g*g,BB/2, where g* is the Lande 5 -factor of 
electron in the QD, gs is the Bohr magneton. The sys¬ 
tem requires spin up electrons that injected from the left 
lead into the dot 1. We take energy levels that satisfy 
£ t > Ml and £4- — £1 Ml, where Ml is inter-dot cou¬ 
pling strength. This energy structure forms a spin condi¬ 
tioned repulsive potential which ensures that the spin up 
electron is forbidden to enter the dot 2 until its spin is 
flipped to be upside down due to its coupling to the mag¬ 
netic atom. As soon as the electron spin changes to be 
down, it would be allowed to pass through the dot 2 and 
collected in the right lead. At the same time, the spin 
of magnetic atom would change from (M z ) to (M z ) + H. 
If a spin down electron is injected from the left lead, it 
directly transports through the double dots without any 
change in the spin of the magnetic atom. However, if 
the left lead is a fully up polarized electron source, each 
passed electron would be connected to the spin change 


of the magnetic atom. In this case, due to spin conserva¬ 
tion, if spin of the magnetic atom is changed from (M~) 
to (. M z ) + n\h, then the number of electrons detected in 
the right lead should equal to and their spins are down 
polarized. 

The magnetic atom can also be driven reversely, which 
is presented in Figs. Qle)-(h). Since the magnetic field is 
turned to be in z direction here, the spin filter only allows 
spin up electrons tunnel through the dot 2. The input 
electrons are required to be in spin down state and out¬ 
put electrons are expected to be in spin up state. Each 
transported electron contributes to the spin of the mag¬ 
netic atom with the value —fi. It yields the impurity 
magnetization orient in the reversal way, which means 
spin of the magnetic atom would be changed from (M~) 
to ( M z ) — n^h with collection of spin up electrons. 

To give a quantitative description to the model, we use 
the following Hamiltonian^^ 


H — Hdl + Hd2 + Hdl2 + Hlead + H tun , (1) 

where the Hamiltonian for dot 1 is 

H d i = £ini + Uihifhii - j e M ■ Si, (2) 

and the Hamiltonian for dot 2 is 

H d 2 = £ 2^2 + U 2 fi2^n 2 i + g*gBB ■ S 2 - ( 3 ) 

The inter-dot tunneling Hamiltonian reads 

Hdi 2 = Kl{fii 2 + n\ 2 ) + Uhifi 2 ■ (4) 


Here, 74 = + h i]r , h it = £^£ 4 , fiii = £4 £ 4 , ™i 2 = 

« 12 t + ^124., « 12 t = C^C2f, hi2\, = C f u C24. Cif (£4) is 
the annihilation operator of spin up (down) electron in 
dot i (i = 1 , 2 ). Si and S 2 represent electron spin in 
the dot 1 and the dot 2, respectively. U is the inter¬ 
dot Coulomb potential. The exchange coupling strength 
between the electron in dot 1 and the magnetic atom 
is given by j e = J|'!/ ; o(’'m )| 2 with exchange integral J 
and the electron ground state wave function ipo at the 
magnetic impurity position j'm ■ 

The left and right electronic leads are described by the 
free electron baths with the Hamiltonian 

Blead — ^ ^ €akC a k a Cak(7 ■ (b) 

k,a]a=L,R 

The dot 1 is coupled to the left leads and the dot 2 is 
coupled to the right lead by 

Htun = ^L £ Lfco- £ l(T + yRCR ka C2a + H.C., (6) 

k,cr k,a 

with the left and right tunneling amplitudes Vj J and Vr, 
respectively. 
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Time evolution of electron transport through the dou¬ 
ble QD system is described by a quantum master equa¬ 
tion which is derived based on the Hamiltonian m and 
the Liouville-von Neumann equation in the Born-Markov 
approximation. Since the transport in our system is a 
process of single electron sequential tunneling through 
QDs and works in the weak tunneling regime, the mas¬ 
ter equation is an effective approach to describe our 
model . 28-30 Equation of motion is given in terms of the 
reduced density matrix p of the system, 

I 


Ap = \ T lUlAAc,PC ler - C la c\aP) 

a 


d 1 ^ 

-q^P = -r^[Hdi + Hd2 + Hdi2, p\ + C-lP + CrP- (7) 

The Liouville super-operators, Cl and C B , acting on 
the density matrix p describe tunneling on the left and 
right side of the double dots, respectively. They are writ¬ 
ten as 


(1 - fL,a)AaPc\ a - c\ a C la p) + H.C .), ( 8 ) 


and 


C-RP = -^'^TrUrAcIoP^o - C2acl a p) + (1 - fR,*)(c2crPcL ~ P) + H.C.), 


(9) 


where the Fermi distribution function in the left lead is 
/l,g = (ea;p[(ei + Uifiig - pL,)/k B T\ + 1) _1 and in the 
right lead is = (exp[(e a + U 2 fi 2s - PFt)/k B T] + 

1 ) _1 , depending on the charging energy Ui conditioned 
by the occupation, of electron with spin 

a. Here, the spin index a is defined that if a =f ( 4 ,), 
then a = 4 , (/). The tunneling rates are spin dependent 
given by r£, = r„(l + P a )/ 2 and r£ = r a (l - P a )/2 
with the current polarization P a = (7^ — )/(l£ + A)- 

Here, 1° is the current with pure spin a on the side of a 
(■a = L,R). The bare tunneling rates can be expressed 
as Tq, = 27rft|f a | 2 N a (e) with the density of states N a (e) 
of electrons at energy e. 

III. MANIPULATION OF A MAGNETIC 
ATOM WITH SPIN M=5/2 

The Hilbert space of the double QD system is gen¬ 
erated by the basic vectors |i)i|m)i|j) 2 , where |i)i and 
|j >2 represent electronic state in dot 1 and dot 2, re¬ 
spectively. Here, i,j = 0 denote the empty state, 
i. j =4,, f indicate occupation states of single electron 
with spin down and spin up, respectively. |to)i is eigen¬ 
state of the impurity spin operator M z with eigenvalues 
?n = — M, —to + 1,..., M. 

First, we consider a typical single magnetic atom 
which displays a spin of M = 5/2 in the QD, such 
as Fe or Mn. These magnetic atoms have six quan¬ 
tized spin states | m), with corresponding eigenvalues 
to = -5/2,-3/2,—1/2,1/2,3/2,5/24 In FigUJa), initial 
state |V , m( 0)) of the magnetic atom is set to be any of 
the six quantized states. In all cases the magnetic atom is 
driven to the final state 15/2). In average value of the Mn 
spin (M z ) = Tr[M z p\ trace is taken over the basic vectors 



FIG. 2: (Color on line) Time evolution of the single impurity 
magnetization as a function of time, (a) The applied external 
magnetic field is B = (0,0, —B). (b) The applied external 
field is B = (0,0, B). The Parameters are pl = 75T, p B = 
-75T, Ei = 0, £2 = 62.5T, g*p B B = 135T, j e = 3T, k B T = 
12.5T, r L = r R = r, p r = o ,m = 5r, u = tor. 


|*)i|m)i|j) 2 . When the impurity spin changes to |5/2), it 
becomes parallel to the injected electron spin and there 
will be no spin flip in the later time. Then the electron 
tunneling should be switched off and the spin down polar¬ 
ized current decreases to be zero as plotted in Fig.[3](b). 
The figures also show that orientation time of the Mn 
spin is a few tens of ns when one takes a tunneling char¬ 
acteristic time of T -1 ~ Ins. It is comparable to the time 
scale of a previously reported optical control^ Since each 
electron contributes to the spin change of the magnetic 
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FIG. 3: (Color on line) (a) The number of spin down electrons 
collected in the right side of the system, (b) Current for spin 
up and down electron s as a function of time. The parameters 
are the same as that in Fig. [2] 


atom with momentum h , the magnetic atom initialized 
in the states |V>m(0)) = | - 5/2), | - 3/2), | - 1 / 2 ), |l/ 2 ), 
|3/2), |5/2) leads to finite electrons collected in the right 
lead with definite numbers ( 7 - 4 ) = 5,4, 3,2,1,0, respec¬ 
tively. It can be seen in Fig.[3la) for corresponding initial 
spin states. The collected electron number is calculated 
using the formula 

(n a ) = - [ Ia(t)dt , ( 10 ) 

e Jo 

where I a {t) is current on the right side of dot 2. The 
current is derived from the charge fluctuating in the two 
dots d((ni) + (ri 2 ))/dt = ~ F^)/e with the 

replacement I<j(t) = JJ). The reversal time evolution of 
the magnetization is shown in Fig. [2] (b). 

IV. MANIPULATION OF A MAGNETIC 
ATOM WITH SPIN M=l/2 

Now, we consider a magnetic ion with spin M = 1/2, 
such as Cu 2+ . This kind of impurity has particular mean¬ 
ing that its maximum magnetization difference is H , from 
—h /2 to H /2 or by inverse. In this case, only a sin¬ 
gle electron is involved in the tunneling. As shown in 
Fig|U)a), when we set initial state of the magnetic atom 
to be | — 1/2), it is transformed into 1 1 / 2 ). At the same 
time, the up polarized current is changed to be down po¬ 
larized current after interacting with the magnetic atom. 
The down polarized current is expected to be collected 
in the right side, which is characterized by a singe elec¬ 
tron with spin down as shown in FigUDb). To show the 
collected charge is really limited, average current as a 
function of time is plotted in the inset of Fig|4jb). The 
down polarized current sharply increases and then disap¬ 
pears slowly with a small fluctuation. The up polarized 
current is negligible weak. From the point of application, 




FIG. 4: (Color on line) (a) Probabilities of the magnetic atom 
in its spin states | — 1/2) and 11/2). The magnetic atom is 
initially in the state | — 1/2). (b) The number of collected 
electrons with spin up and spin down corresponding to the 
case (a). The parameters are pz, = 70T, fiR = —70F, ei = 0, 
£2 = 57.5F, sVbB = 125F, j e = 2F, k B T = 10F, V L = V R = 
T, P L = 1, P R = 0, Ml = 3F ,U= 10F. 


it is a very good system for date storage, where only sin¬ 
gle electron of current and the two states of the magnetic 
atom is correlated. 

The small fluctuation in the current is occurred due to 
the back action from the system. Since charge transfer 
through the double dots requires electron spin flip, the 
amplitude of current is proportional to rate of electron 
spin change. The rate of spin change is determined by 
the strength of spin coupling between the magnetic atom 
and electrons. Therefore, we can deduce that the oscilla¬ 
tions observed in Figj4] is contributed from the coherent 
coupling between the magnetic atom and an individual 
electron. Indeed, character of the small oscillation can 
be tuned by changing the spin coupling strength j e . Ac¬ 
tually, the back action effect can also be seen in the for¬ 
mer situation for M = 5/2. However, considering six 
spin states of the magnetic atom are involved in the ex¬ 
change interaction and the larger couplings j e and Q are 
taken, the small oscillation during the evolution of the 
spin states and current is too smooth to be observed. 

The above results are described in a relatively ideal 
situation. At the end of this section, let us talk about 
some more practical cases. Fig.[5])a) reveals that increase 
of the inter-dot tunneling strength excites a spin up elec¬ 
tron from the dot 1 into the excited level 4 of the dot 2 . 
It leads to the unexpected spin up electron leakage from 
the double dot system into the right lead. To restrain the 
electron leakage through the excited level e-j-, we suggest 
that one can take a relatively small inter-dot coupling hfl 
which is required to be much smaller than the decoupling 

£ t - £1. 

The system is also sensitive to polarization of the right 
lead. In Fig. [5jb) , it is illustrated that when the polar¬ 
ization is not pure, the number of electrons collected in 
the right lead is larger than the expected value. Since 
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FIG. 5: (Color on line) (a) Number of spin up electrons col¬ 
lected in the right side of the double QD with different inter¬ 
dot coupling strength, ksT = 10r, Pl = 1. (b) Number of 
spin down electrons for different spin polarization in the left 
lead, ksT = 10F, hQ. = 3F. (c) Number of collected spin 
down electrons at different temperature, MI = 3r, P L = 1. 
The rest parameters are the same as that in Fig. [4] 

non-purely polarized current contains electrons with dif¬ 
ferent spin orientations, thereby some electrons could not 
be blocked by dot 2. 

The negative number of electrons in Fig. [5jc) implies 
that electrons in the right lead have certain probability 
to flow into the QDs at the beginning of tunneling due to 
the thermal excitation, since the QDs is empty initially. 
After the system reaches steady state, the higher the tem¬ 
perature, the larger electron distribution would be in the 
QDs. The electron staying in the QDs comes from the 
leads. Therefore, net number of electrons collected in the 
right lead decreases when temperature increases. 

V. DISCUSSIONS 

In experiment, to realize electron transport through a 
magnetic atom doped II-VI semiconductor QD is still a 
challenge. There are some well developed experimental 
backgrounds, for instance, the spin dependent electron 
transport through a quantum well that containing dilute 
magnetic material^ and resonant tunneling through a 
CdSe self-assembled QD with Mn ions— In addition, an 
attempt to orient a Mn spin using the charge that trans¬ 
ported from a neighboring QD has been successfully re¬ 
alized^ 

To create a local magnetic field is a key technical prob¬ 
lem for the experimental implementation of our proposal. 
The left lead is a polarized electron source, so an exter¬ 
nal magnetic field may be required on it. In addition, 
the right QD needs a magnetic field to create a Zeeman 
splitting for an electron. At the same time, these mag¬ 
netic fields must have negligible affect to the left QD and 


the right lead. Therefore, some local magnetic fields are 
necessary in the nanostructure. To this end, a proper 
quantity of external magnetic field can be exerted on 
the left lead for the generation of spin polarized electron 
source. As the QD is a zero dimensional nanoscale sys¬ 
tem, our suggestion is using a magnetic grain, such as Co 
grain) 33 ! 34 to create magnetic field on the right QD. To 
explore its effectiveness we estimate the field strength of a 
typical Co grain which consists of hundreds of atoms. We 
assume the Co grain is a uniformly magnetized sphere. 
Then, the field of the grain can be obtain from the for¬ 
mula B = |^oM, /r 0 is permeability of free space and M is 
magnetization of the grain. 35 We take that average mag¬ 
netic moment per Co atom is 1.5/j.b- In real materials the 
local magnetic moment is always larger than this value. 33 
We take diameter of a Co atom is about 0.15nm. Then, 
magnetic field of the Co grain is estimated to be not less 
than B = 6.6 T. By increasing size of the Co grain, more 
strong field can be obtained. In a CdTe semiconductor 
QD, the corresponding Zeeman splitting of an electron 
reaches \g* llbB\ = 0.64meV, where the (/-factor is given 
by g* = —1.67 in this material. 3 S To guarantee the left 
QD and the right lead is not effectively influenced by the 
surrounding magnetic field, magnetic shielding may be 
applied here to protect the spreading field. One kind of 
the magnetic shielding materials is superconducting chip 
which expels magnetic field via the Meissner effect. An¬ 
other kind shielding material is certain high permeability 
metal alloy which, in contrast the superconductor, draws 
the held into themselves. Taking an example with the 
high permeability shielding with permeability /r, held in 
the shielded volume is 2 {b 3 ~^a 3 )^ H mes smaller than the 
outside held.— For the convenience of quantitative esti¬ 
mation, the shielding material here is assumed to be a 
spherical shell with inner radius a and outer radius b. 
When the permeability /r is large enough, good shield¬ 
ing of the held in the shielded area can be achieved. In 
practice, as the QD couples to an electronic reservoir and 
another dot, the shielding material may be not absolutely 
closed with a lower efficiency. 

In the Hamiltonian of our model, spin exchange inter¬ 
action between the two QDs is not considered. In fact, 
when each quantum dot contains one electron, the ex¬ 
change interaction does not play important role for the 
whole system. The reason is that the electron in dot 2 
has dehnite spin direction. It has to gain large energy 
to change spin. However, it is hard for the electron in 
dot 1 to provide the large energy. Besides, the system is 
robust against the double occupancy in dot 1. Since cur¬ 
rent injected from the left lead is assumed to be pure spin 
polarized, two electrons with different spin states in dot 
1 does not break the spin conservation as soon as dot 2 
guarantees to output electrons with pure spin. Whereas, 
double occupancy in dot 2 induces spin leakage and, as 
a result, the correlation between electrons and the mag¬ 
netic atom becomes weak. Even though, considering the 
high Coulomb blockade effect in either of the QDs, double 
occupancy in any dot is negligible small in our model. 
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As an information processing system, its characteristic 
times are very significant. There are two kinds of critical 
times, one of them is manipulation time r s of the system, 
another is lowest life-time bound t 7 ; of these information 
carriers, such as the magnetic atom and electrons. The 
manipulation time of the system indicates a time range 
during which an electron is emitted from the left lead and 
then collected in the right lead, at the same time control 
of the magnetic atom is accomplished. It is clear that 
spin life times of the information carriers are required to 
be, at least, longer than the manipulation time of the 
system, i.e. t, > t s . As mentioned in the introduction, 
a single Mn spin relaxation time from 1 [is to 0.4 ms 
in CdTe QEUJI and electron life time from 50 ns to 1 
[is in II-VI semiconductor systems are reported. Even 
in a QD including certain charges, the Mn atom relaxes 
in a time scale of about 100 nsM As shown in section 
III and IV, manipulation time of the system is about 
50 ns for a characteristic time of the electron tunneling, 
T _1 ~ Ins. Furthermore, there is still a space for re¬ 
ducing the manipulation time of the system. Actually, 
properly increasing the dot-reservoir coupling (electron 
tunneling rate) or inter-dot coupling would improve the 
rate of control process. In this case, a more lower bound 
of the spin life time is allowed. 


VI. CONCLUSIONS 

In the system, change of magnetic impurity spin is cor¬ 
related to spin state and number of single electrons that 
tunneling through the two QDs. Based on this principle, 
we give the following predictions: (i) Our model works 
as an electron source in which number of emitted elec¬ 
trons can be determined beforehand by setting an initial 
state of the magnetic atom or using a magnetic atom with 
certain spin M. In particular, a single electron emitter 
is available using a magnetic atom with spin 1/2. (ii) 
The number of polarized transferring electrons can be 
recorded in the spin state of the magnetic atoms, (iii) 
The change in spin state of the magnetic atom should 
be detected by counting the number of electrons that 
emitted from the double QDs. (iv) Spin state of the 
magnetic impurity can be controlled, in principle, by in¬ 
jecting suitable number of spin polarized electrons, and 
this controlling is reversal. 
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